Solar induced growth of silver nanocrystals 
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The effect of solar irradiation on plasmonic silver nanocrystals has been investigated using Trans- 
mission Electron Microscopy and size distribution analysis, in the context of solar cell applications 
for light harvesting. Starting from an initial collection of spherical nanocrystals on a carbon film 
whose sizes are log-normally distributed, solar irradiation causes the nanocrystals to grow, with 
one particle reaching a diameter of 638 nm after four hours of irradiation. In addition some of the 
larger particles lose their spherical shape. The average nanocrystal diameter was found to grow 
as predicted by the Ostwald ripening model, taking into account the range of area fractions of the 
samples. The size distribution stays approximately log-normal and does not reach one of the steady- 
state size distributions predicted by the Ostwald ripening model. This might be explained by the 
system being in a transient state. 

PACS numbers: 



I. INTRODUCTION 

Metallic nanocrystals are desirable for their use in so- 
lar cell applications, where they can be used to increase 
efficiency by means of the localized surface plasmon reso- 
nance. This collective oscillation of the valence electrons 
results in a local electromagnetic field at the nanocrys- 
tal surface and in wavelength selective photon absorption 
and scattering. By dispersing metal nanocrystals on the 
surface of an optically thin solar cell, the localized sur- 
face plasmon resonance will scatter the light further into 
the solar cell, thereby increasing the absorption and thus 
the efficiency. 

Silver (Ag) nanocrystals are particularly popular, be- 
ing relatively cheap and easy to produce by a variety 
of processes. In addition, they do not need to be de- 
posited and heat-treated, which can be damaging for 
the solar cells. Using Ag nanocrystals, PillaP and Beck 
et al.'^ demonstrated a significant increase in absorption 
in wafer-based silicon solar cells. However, nanocrystals 
used in solar cell applications will be exposed to elevated 
temperatures, which can make them unstable. A solid 
understanding of the underlying growth and coalescence 
mechanisms is therefore of vital importance. 

The Ag nanocrystals in this work have been made 
by a chemical synthesis of metallic nanocrystals, using 
chemical reduction of silver salts by sodium borohydride. 
This has shown to be a simple, economical, and pop- 
ular method that can be applied in large scale as re- 
quired for industrial application^^. The surface plas- 
mon resonance depends on the shape of the nanocrystal 
surface, size, spatial arrangement, and configuration of 
the nanocrystal^. It is therefore very important to in- 
vestigate these features in detail. We have previously 
studied the particle size distribution, structure^ and op- 
tical properties of colloidal Ag nanocrystal^^. These 



nanocrystals, however, must remain stable during solar 
irradiation and heat. In this paper we study the effect 
of solar irradiation on Ag nanocrystals on a carbon film, 
and investigate the growth rate and particle distribution 
with solar irradiation time. 



II. EXPERIMENTAL 

The Ag nanocrystals were made by a wet chemical reduc- 
tion synthesis mixing a silver nitrate solution (AgNOa, 
99.9% purity, Qualigens Fine Chemicals) with a highly 
reducing solution of sodium borohydride (NaBH4, 95% 
purity, Merck) without further purification. All equip- 
ment was thoroughly cleaned by soaking in ethanol and 
washed with distilled water. A 0.6 mM AgNOs precur- 
sor solution was prepared using deionized II2O, together 
with a 1.2 mM aqueous NaBIl4 solution. In order to 
prevent agglomeration, the solutions were kept ice cold 
and the mixture was stirred continuously. Ag nanocrys- 
tals were formed by adding the aqueous AgNOa drop by 
drop to the highly reducing solution of NaBH4. Varia- 
tions in particle size and distribution were acquired by 
varying the amount of AgNOa (aq) relative to NaBH4 
(aq) in the final mixture. The samples studied in this 
work have a AgN03/NaBH4 molar ratio of 2/25. 

Six Transmission Electron Microscopy (TEM) samples 
were prepared by adding a drop of the final solution on a 
holy carbon film supported by a Cu grid. In order to ob- 
serve light induced efi^ects, five of the TEM samples were 
placed under a Xenon lamp, irradiating the wafers with 
an intensity of 100 mW/cm^ (1 Sun). While the samples 
started off at room temperature, the temperature slowly 
increased from 55.7°C after 15 min. to 64.3°C after 240 
min. The lamp was placed 35 cm from the TEM sample 
grids. The remaining sample, which we shall refer to as 
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sample 1, was kept aside. One by one the other samples 
are then taken away from the Xenon lamp, sample 2 af- 
ter 15 min., sample 3 after 30 min., sample 4 after 60 
min., sample 5 after 120 min., and finally sample 6 after 
240 min. The Ag nanocrystals were then observed by 
High-Resolution Transmission Electron Microscopy in a 
200 keV JEOL 2010F microscope with a Gatan imaging 
filter and detector. 



III. THEORETICAL 

A. Ostwald ripening 

Ostwald ripening is a thermodynamically-driven spon- 
taneous process, in which small particles dissolve and 
through diffusion redeposit onto larger particles in an 

effort to decrease the total free surface energy. 

While first qualitatively described by OstwalcP'^, the 
theory made a major leap forward when Lifshitz and 
Slyozo\ffi' and WagneiQ^ independently derived quanti- 
tative predictions about the distribution and evolution 
of the particle sizes. The LSW model, as it is now called, 
was derived under several assumptions, including: 

1. The particles are spherical. 

2. The distances between the particles are infinitely 
large compared to the size of the particles. 

In practice, however, particles are not always spherical, 
and infinite distances are physically unrealistic; the dis- 
crepancy is characterized in terms of the volume fraction 
(f> (or area fraction in 2D), which is defined as the volu- 
metric density of the particles. It is therefore no surprise 
that the size distribution predicted by the LSW model 
corresponds to virtually none of the size distributions 
encountered in experiments, which are generally broader 
and more symmetric!!^. 

To address these shortcomings several extensions to 
positive volume fractions were proposed, in cludin g an- 
alytical calculation^^ and numerical modelJ^^'^, some 
extending to nonspherical particle^^^. For a comprehen- 
sive survey on the progress in Ostwald ripening theories, 
see the paper by BaldmarP. 

One common feature of these extensions is the pre- 
diction that the average particle radius (i?) follows the 
growth law 

{Rf = {Rof + k^t, 

where (i?o) is the average particle radius at i = 0, and 
the kinetic coefficient kip depends on the area fraction (p 
of the sample. Although k^ almost always increases with 
(/>, there are known situations where it decreases with 
insteacff^. 

Since some of our particles become moderately aspher- 
ical over time, they do not have a well-defined radius. 
We simplify our analysis by assuming the particles are 



spherical, and then computing their diameter d from 
the measured cross-sectional area A from the formula 
A = 7r((i/2)^. As the diameter of a spherical particle 
is twice its radius, we expect the average diameter (d) to 
satisfy a growth law 

= (do>'+8M- (1) 
Numerical simulations by Fan et suggest that 

fco 25 ~ 0.810, fco 5 ~ 2.548, fco 75 « 8.207, ko 9 w 24.487. 

(2) 

B. The log-normal distribution 

In our previous work on the growth of colloidal Ag 
nanocrystals, the particle sizes were shown to be log- 
normally distributecP. The discussion gave insight into 
the growth and structure of the formation process and 
provided a close fit with the experimental observations 
of the particle sizes. 

The log-normal distribution with parameters /i and a 
has probability density function 

1 (In x-m)^ 

p{x]fi,a) = j=e 2.2 ^ a; > 0. (3) 

xay/ Ztt 

The parameters /i and a coincide with the mean and stan- 
dard deviation of the natural logarithm of the random 
variable X, which is by definition normally distributed. 
The random variable X itself has mean 

m = e'^+'^'Z^ (4) 

and standard deviation 

= e^+"'/Ve^' - 1. (5) 

C. Parameter estimation 

In our previous paper we chose to smoothen the overly 
erratic particle size densities by computing moving 
averages^. Since there is no clear choice for the size of the 
sampling window for the data in this paper, we instead 
choose to group the diameters into bins of each 1 nm in 
width. 

The probability that the random variable X assumes 
a value between d — 1/2 and d -I- 1/2 is 

l'd+1/2 

Pd;t,,cr ■= / p{x;fi,a)dx. 

Jd-1/2 

Suppose now that sample i comprises Ni particles of 
which Ni,d are of size between d— 1/2 and d-l- 1/2. Then 
the log-likelihood i^ 

^Ar,.rflog(Pd;^,,). (6) 
d 
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The maximum likelihood estimators /i, of /i, cr are found 
by maximizing this expression numerically. 

Plotting, for each sample, the mean nanocrystal diame- 
ter against the irradiation time suggests that there exists 
a relation ([T]) between these two quantities. For a set 
of uncorrelated data with varying uncertainties, as is the 
case here, it is unreasonable to assume that every data 
point should be treated equally in estimating an underly- 
ing relation. Instead, a best linear unbiased estimator is 
obtained by weighing each data point by the reciprocal of 
its variance^^. This weighted least square method is well 
known, and was briefly described in a similar context in 
our previous worlP. 



IV. RESULTS 

The size and shape of the particles has been studied with 
TEM, before and after solar irradiation. Initially the par- 
ticles are spherical and small, with an average diameter 
of ^6 nm. During the solar irradiation most of the par- 
ticles grow, at the cost of the smaller particles. Some of 
the particles grow very large, the largest observed particle 
having a diameter of 638 nm after four hours of solar irra- 
diation, and lose their spherical shape. Not all particles 
grow, however, as the sample irradiated for four hours 
still contains nanocrystals with the minimal observable 
diameter. Electron Energy Loss Spectroscopy revealed 
very little oxidation of the Ag nanocrystals. Only the 
outermost surface of the largest nanocrystals seems to 
be affected. 



A. Particle distribution 

For each sample at least six TEM images were taken 
and analyzed; see Figure l^for one image of each sam- 
ple. The program ImageJ^S was used to compute the 
area fractions in the samples, and to accurately count 
the particles and estimate their diameters. The average 
area fraction in each sample varied from 0.240 to 0.463, 
while the number of counted nanocrystals varied from 
257 to 1499. Since the bigger nanocrystals in sample 6 
(partially) overlap many of the smaller nanocrystals, the 
nanocrystals in these TEM images were separated man- 
ually before applying Image J. A qualitative check was 
carried out by manually measuring the diameters for a 
smaller number of representative nanocrystals. Figure [2] 
shows, for each sample, a histogram of the diameters, 
the bins being of 1 nm in width and centred at integral 
diameters. 

Next, the maximum likelihood estimators /I, ct are 
found by maximizing Q numerically, and the mean and 
standard deviation are computed from Q and ([s]) . These 
values are listed in Table |Tj and the estimated log-normal 
densities are plotted in Figure [2] 




FIG. 1: TEM images of the Ag nanocrystals after solar irra- 
diating for 0, 15, 30, 60, 120, and 240 minutes. Note that the 
images of samples 5 and 6 are taken at a different resolution. 
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TABLE L For each of the samples i = 1, 2, . . . , 6, the table 
lists the irradiation time ti in minutes, the number Ni of 
counted nanocrystals, the parameters /ii and ai of the fitted 
log-normal distributions, and the corresponding expectation 
value rrii and standard deviation sdi. 



B. Particle growth 

Table |l] suggests that the mean nanocrystal diameters 
rrii grow with the irradiation time ti. For a quantita- 
tive analysis we plot these quantities against each other 
in Figure |3j with error bounds the corresponding stan- 
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FIG. 2: The figure shows the relative number of Ag nanocrys- 
tal for each diameter, after solar irradiation for 0, 15, 30, 60, 
120, and 240 minutes. Also drawn are fitted log-normal densi- 
ties, whose parameters fi and a were estimated by maximizing 
the log-likelihood (|6| numerically. 




60 120 180 

Irradiation time (min.) 



FIG. 3: Average diameters against the irradiation times, with 
the error bar extending one standard deviation in each direc- 
tion. The dashed curve is the growth law ([T| estimated by 
the weighted least square method. Drawn dotted are the pre- 
dicted curves corresponding to the minimal and maximal area 
fraction. 



V. DISCUSSION 

A. Particle distribution 

As shown in a previous worlP^ the log-normal distribu- 
tion provides a close fit for the diameters of the nanocrys- 
tals in sample 1. After irradiating the samples for various 
times, the log-normal distribution continues to be a good 
fit for the diameter distributions for samples 2-5. In 
sample 6 the log-normal distribution no longer provides a 
good fit, and the steady-state size distributions predicted 
by the Ostwald ripening model do not fit either. How- 
ever, the peak is becoming more pronounced, and after 
scaling each empirical density by the average diameter, 
the tail seems to be moving inwards. These features are 
consistent with the system being in a transient statcP. 

The evolution of the particle size distribution under 
the Ostwald ripening is notoriously slow, and the time 
required to reach a steady state generally depends on 
the initial distributiorP3! and the area fraction^'^. Com- 
pared to sample 1, the average diameter in sample 6 has 
changed by a factor of mg/mi = 12.572/6.134 ^ 2.05. 
This is typically too little to reach the steady stat^^Sl^ 
another indication that the observed distribution is in 
a transient regime. More data, in particular data for 
longer irradiation times, are needed to confirm this phe- 
nomenon. 



dard deviations sdi. Drawn dotted are the curves ([!]) 
corresponding to the minimal area fraction (w 0.25) and 
maximal area fraction (w 0.5), with (do) ~ 6.152 the 
empirical mean diameter of sample 1 and k^j, as in |^ 
predicted by the numerical simulation by Fan et al!^. 
Applying the weighted least square method, with weights 
the reciprocals of the variances sdf, yields a fit of ([l]) with 
(do) = 6.429 and fc^ — 0.881. The corresponding curve 
is plotted (dashed) in the same figure. 



B. Particle growth 

With longer solar irradiation the particles grow and be- 
come aspherical. Such aspherical particles are common 
for higher area fractions^. The growth of the average 
diameter of the nanocrystals fits the growth law (I]) pre- 
dicted by the Ostwald ripening model. Because the area 
fraction varies from sample to sample, the kinetic coef- 
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ficient takes on different values as well. As can be seen 
in Figure [3j however, the least square estimate fits well 
between the corresponding range of kinetic coefficients 
predicted by Fan et all^^. 

As seen in our previous worl5P, the increase in particle 
size will yield an increase in the wavelength of maximum 
scattering. This will increase the efhciency of the device 
initially, since the nanocrystals investigated in that pa- 
per have a maximum scattering peak at relatively low 
wavelengths (e.g. ~ 480 nm in silicon nitride), whereas 
the optimal efficiency for the solar spectrum is attained 
above 600 nm, for instance at 700 nm. Since the device 
would be unstable, however, its efficiency would change 
with time and would eventually decrease again. 

VI. CONCLUSION 

The effect of solar irradiation on Ag nanocrystals has 
been investigated by TEM and size distribution analysis. 
Solar irradiation causes the nanocrystals to grow and to 
lose their spherical shape. After four hours of solar irradi- 



ation, the largest observed particle has a diameter of 638 
nm and is still surrounded by small nanocrystals. The 
average particle diameter (c?) can be approximately de- 
scribed by the growth law ([!]) , with initial average diame- 
ter (do) ~ 6.429 and kinetic coefficient « 0.881, where 
d is in nanometers and t is in minutes. While the particle 
size distribution initially stays log-normal, it eventually 
starts to deviate from the log-normal distribution. The 
resulting distribution does not fit the steady-state distri- 
butions predicted by the Ostwald ripening model either, 
which might be explained by the system being in a tran- 
sient state. 
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